Effect of different bending shapes on thermal properties of flexible light-emitting diode filament
Wang Liping1, Li Wenbo3, Xu Yichao1, Yang Bobo2, Shi Mingming2, Zou Jun2, †, Li Yang1, ‡, Qian Xinglu1, Zheng Fei1, Yang Lei3
School of Materials Science and Engineering, Shanghai Institute of Technology, Shanghai 201418, China
School of Science, Shanghai Institute of Technology, Shanghai 201418, China
Zhejiang Emitting Optoelectronic Technology Co., LTD, Jiaxing 310000, China

 

† Corresponding author. E-mail: zoujun@sit.edu.cn liyang123@sit.edu.cn

Project supported by the National Natural Science Foundation of China (Grant No. 51302171), Science and Technology Commission of Shanghai Municipality, China (Grant No. 14500503300), Shanghai Municipal Alliance Program, China (Grant No. Lm201547), Shanghai Cooperative Project, China (Grant No. ShanghaiCXY-2013-61), and Jiashan County Technology Program, China (Grant No. 20141316).

Abstract

Heat dissipation is an important part of light-emitting diode (LED) filament research and has aroused constant concern. In this paper, we studied the thermal performance of flexible LED filament by numerical simulation and through experiment. The heat dissipation characteristics of spring-like structure flexible LED filament were computed by finite volume method, and it was found that the chip junction temperature was closely related to the pitch and the bending radius. The effect of inclination angle of lighting LED filament was discussed because it is relevant to the spring-like structure flexible LED filament in geometry. The results demonstrated that the temperature of the filament increases as the inclination angle improves.

1. Introduction

Light-emitting diode (LED) lighting technology has upgraded and improved so quickly that LEDs are now widely used in a wide variety of applications.[13] Phosphor-converted light-emitting diodes are the most popular white LED products owing to their high efficiency, low cost, and high reliability.[48] In 1997, Nichia made the first commercial white LED all over the world in the manner of light conversion, which was fabricated by coating yellow phosphors on the chip of blue LED. The yellow light emitted by the phosphors was mixed with the blue light from chips into white light.[9] As chip power and luminous efficiency increase every year, LEDs have become a new type of light source. In April 2013, Philips Electronics N.V. announced the release of new products whose luminous efficiency was up to 200 lm/W, which is much higher than the incandescent lamp (100 lm/W) and metal halogen lamp (15 lm/W). Temperature control is one of the most important technologies in modern LEDs.[1013] It was reported by Arik et al.[14,15] that the heat accumulated in chips raised the temperature of PN junction, and this phenomenon led to shorter service life, low color-rendering index, and other issues. Therefore, the heat dissipation problem of LED light sources has become an important factor influencing the LED luminous efficiency and luminous decay.

The packaging methods and materials of LEDs, and the environment that they are used in have a significant influence on the luminous effect. Improving the heat dissipation is a viable approach to enhance the performance of LEDs. There has been a tremendous change in filament technology from tungsten filament to LED filament. During the Guangzhou international lighting exhibition in 2013, only two companies (Shanghai Oppel-LED Co. Ltd and Hangzhou Jingyang Co. Ltd) exhibited LED filament bulbs. In contrast, during the Hong Kong international lighting fair in 2014, there were more than 200 companies exhibiting all kinds of LED filament bulbs. It was reported that more than 600 enterprises produced the LED filaments. The filaments adopted high voltage technology (HVLED) and small current driving technology, which would reduce the heating degree of the light source greatly. Most studies about filament heat dissipation have only focused on the rigid filament, the substrate of which is too hard to bend casually. At the same time, the thermal properties concerning the horizontal placement angle have not been discussed in detail.

Computer numerical simulation technology is an important method for thermal simulation and integrated circuit design.[16] According to Christensen,[17] the heat dissipation performance of LED modules that include a single Lambertian LED and a single chip-on-board (COB) packaged ceramic substrate under different heat dissipation conditions was investigated using the finite element method. In this study, a computational model was used to simulate the Philips (E27ES) LED filament bulb. However, most studies on the bulb temperature have only focused on the filling gases and filament substrates, and no single study about the thermal properties on filament shape has been researched.

In this paper, we will analyze the mathematical model of heat conduction and the flexible filament model of convection computer through numerical simulation method, and also the equivalent thermal resistance. The influence of the filament pitch on heat dissipation is also analyzed. At the same time, the temperature of the filament is tested with an infrared temperature measuring instrument, and the result is compared with the numerical simulation.

2. Mathematical modeling
2.1. Experimental setup

The manufacturing process flow of flexible LED filament with a length of 100 mm is shown in Fig. 1. An average of 100 flip-chips was distributed on the copper substrate, the size of which was 200 μm × 510 μm × 250 μm. SY-6021 silicone A (Shenzhen Shenghuayang Electronic Material Co. Ltd), SY-6021 silicone B (Shenzhen Shenghuayang Electronic Material Co. Ltd), Lu3Al5O12: Ce phosphor (particle size 13 μm, emission band 530 nm, Yinghe in China), (Sr, Ca)AlSiN3: Eu phosphor (emission band 640 nm, Yinghe in China), and (Sr, Ca)AlSiN3: Eu phosphor (emission band 628 nm, Yinghe in China) were used in this study. The insulation was located in the middle of the paste to prevent short circuits of the electrode. When these procedures were finished, a flexible filament could be obtained. In Table 1, the component dimensions of the multiple chip modules and the thermal properties of the encapsulated materials are shown.

Fig. 1. (color online) Graphic of the process flow of flexible LED filament.

Flexible filaments without and with injection currents are shown in Figs. 2(a) and 2(b), respectively. A single chip was extracted to the model according to the symmetry of the package structure, and the temperature field and heat transfer process were calculated when the flexible filament is working. We assume that 80% of the rated power of the filament input is converted into heat energy.[47]

Fig. 2. (color online) Flexible filaments (a) without and (b) with injection currents.
Table 1.

Filament module component size and thermal property parameters of packaging materials.

.
2.2. Geometry

Insulated-metal-substrate (IMS) was used as the power electronic substrate of this LED filament. The flip chips were attached on the substrate through Au–Si eutectic bonding. Finally, the flip chips were attached on the substrate through solder and protected by the phosphor glue. A glass bar was used to support the shape of the flexible filament.

2.3. Model establishment

The following conditions were fulfilled to simplify the calculation process for numerical analysis:

The calculation has a major emphasis on the steady state of thermal system.

Air density is calculated by the approximate method of Boussinesq.[18]

Owing to the negligible difference when using temperature-dependent thermal conductivity, the thermal conductivities of all materials are independent of temperature.[19]

All of the physical properties of air except density are regarded as fixed values.

Conservation of mass, energy, and momentum is the basic equation to solve the energy equation describing the heat and mass transfer of the system. In the calculating model, the mass conservation equation is

where ρ is the air density, and u, v, and w are the speeds in x, y, and z directions, respectively.

The momentum equation is as follows (−y is the direction of gravity):

where P is the pressure of gas, μ is the viscosity of air, g is the acceleration of free fall, and β is the air expansion coefficient. The energy conservation equation is
where k is the thermal conductivity, cp is the specific heat at constant pressure.

The air density is calculated by the approximate method of Boussinesq

where ρ0 is the air density at room temperature.

Heat can be delivered by convection, radiation, and thermal conduction. Heat is delivered by thermal conduction mainly in solid state. In the interfacial region, heat is delivered in the way of convection. In fluids, both thermal conduction and convection exist at the same time. Radiation could be emitted from solid surface to the colder outside by means of electromagnetic wave. The packaging model and boundary conditions are shown in Fig. 3. The other faces are open boundaries for environmental temperature and environmental pressure. The junction temperature of chips is taken as the monitoring point. The other coefficients of heat transfer could be calculated under the coupling condition between solid and fluid. Assuming that Ta is the environmental temperature at 25 °C, τ is the direction of heat transfer, then the Fourier's law is

where k is the thermal conductivity, and A is the surface area of the filament.

Fig. 3. (color online) Flexible filament module and boundary conditions.

Convection exists in the region between filament model and air. Newton cooling formula is

where h is the convection heat transfer coefficient.

2.4. Numerical procedure and validation

Figure 4 shows the numerical model of solder, single chip, substrate, and phosphor layer of the package module. These models are designed in real size. The temperature field and heat transfer process at working state are calculated.

Fig. 4. (a) Solder, (b) chip, (c) substrate, and (d) phosphor layer.
2.5. Mesh and boundary conditions

As shown in Fig. 4, there are four sub-domains included in the computational domain; i.e., chips, aluminum substrate (Al substrate), phosphor layer (glue), and also gas outside the filament. The outside cuboid enveloping the filament is included to simulate the heat change between the filament and the ambient. Heat conduction, natural convection, and radiation heat transfer are involved in the problem. Standard heat conduction equation is used for heat transfer in solid domains such as chips, substrate, and glue. Navier–Stokes equations with buoyancy force term are used for natural convection outside the filament. Figures 5(a) and 5(b) show the multi-block unstructured meshes which are generated by FLOEFD 15.0 in all domains. Chips are distributed on the filament substrate evenly. Figure 5(c) shows the physical model of the computational domain clearly according to the side view.

Fig. 5. (color online) Numerical simulation model: (a) mesh, (b) side view, and (c) the physical model of the computational domain.
3. Results and discussion
3.1. The effect of a spring-like filament pitch

Figure 6 shows the thermal distribution of the spring-like LED filaments at different pitches. In Fig. 6, we have demonstrated that the highest filament temperature decreases widely by increasing the stretching pitch while the lowest temperature is almost constant. Compared with straight filament, the heat generated by the multi-chips of cylinder screw distribution is more severe owing to the closer distribution, which prevents heat from spreading outward.

Fig. 6. (color online) Thermal distribution of filaments at different stretching pitches by simulation.
3.2. The effect of a spring-like filament bending radius

The result of the thermal distribution of the filament at different stretching pitches by infrared test is shown in Fig. 7. It can be seen that the highest filament temperature decreases widely by increasing the stretching pitch while the lowest temperature is almost constant, which agrees with the simulation result as shown in Fig. 6. The thermal conductivity of the substrate is higher than that of the silica. Thus, it is expected that this will result in the higher temperature of substrate than that of phosphor silica.

Fig. 7. (color online) Thermal distribution of filaments at different stretching pitches by infrared test.

By combining Figs. 6 and 7, we can draw the conclusion that the highest temperature point is the upper part on the filament, but not the top. The larger the pitch, the longer the high temperature part, and the lower the air heat density. This indicates that a larger pitch results in a lower filament temperature.

Figure 8 further illustrates the highest temperature and lowest temperature of the filament by simulation and infrared test under different pitches (measured at a typical injected power of 3.5 W) from 11 mm to 19 mm. It is evident that the highest temperature of both simulation and infrared test decreases with the increase of the filament pitch, while it is opposite for the lowest temperature. This indicates that a larger pitch benefits the heat dissipation, whereas it simultaneously increases the heat exchange area surface between air and filament. Furthermore, the ambient temperature is responsible for the lowest temperature of filament. Some reasonable reasons behind these results are provided.

Fig. 8. (color online) Highest temperature and lowest temperature of filament by simulation and infrared test under different pitches.
3.3. The effect of the filament inclination angle

The spiral line can be understood as a complex movement of the particle, together with circular motion and straight motion, as shown in Fig. 9. The ratio of pitch to radius eventually affects the helix angle. The temperature distribution of simulated straight filaments at different inclination angles is depicted in Fig. 10. It can be seen that the highest filament temperature increases widely by increasing the inclination angle, while the lowest temperature is almost constant.

Fig. 9. (color online) Sketch of spiral filament expansion.
Fig. 10. (color online) (a) Schematic diagram of the inclination angle. (b)–(f) Thermal distribution of filaments at different inclination angles by simulation.

Figure 11 shows the corresponding relation of the highest temperature under different inclination angles, in which the black squares correspond to the maximum temperatures in Figs. 10(b)10(f).

Fig. 11. Corresponding relation of the highest temperature under different inclination angles.
3.4. The thermal resistance of the LED filament model

Thermal resistance is an important parameter to measure the performance of heat dissipation. The model can be calculated under the total power of 3.5 W through Eq. (9). The thermal resistance of the filament is the thermal resistance of multi-chips ranging in parallel. The total thermal resistance is reduced with the increase of the number of chips. Multi-chips in parallel increase the contact area, resulting in the reduction of the thermal resistance. It means that at a constant input power, the more chips, the better the heat dissipation, which can be calculated through Eq. (10).

There are 4 basic steps of the heat transfer: (I) from chips to glue; (II) from chips to substrate; (III) from glue to ambient; (IV) from substrate to ambient. Each step is closely related to one or more thermal resistance models.

There may be thermal contact resistance at the interface between chips and glue, chips and substrate in step (I), and also a spreading thermal resistance when heat is transferred from chips to glue or from chips to substrate. It is difficult to accurately analyze the spreading thermal resistance owing to the complex structure of the filament. The thermal resistances concerning step (I) are not covered in the discussed thermal resistance network and, therefore, the suggested model is only capable of predicting the average filament temperature.

When heat is transferred from the glue or substrate to the ambient in step (II), it is partly transferred by natural convection and partly through radiation, associated with the convective resistance and radiative resistance, respectively. There is only conduction in the heat transfer process of glue and substrate, corresponding to the conductive resistance.

The thermal resistant network, beginning from the junction temperature and ending at the ambient temperature, is shown in Fig. 12, where Ta is the ambient temperature, RNC is the convective resistance, Rrad is the radiative resistance, Rglue, cond is the conductive resistance of glue, Tj is the temperature of PN junction, RAl, cond is the conductive resistance of substrate, and Rc − p is the thermal resistance of chip.

Fig. 12. Schematic diagram of thermal resistance model for LED filament.

The total thermal resistance R is

where δ is the average thickness of substrate and glue.

4. Conclusion

In summary, the shape of flexible LED filament has a significant influence on the heat dissipation of spring-like structure flexible LED filament. The junction temperature of the chip is closely related to the pitch and the bending radius. The results demonstrate that the temperature of the filament increases as the inclination angle improves. It can also be calculated that if more chips are distributed at a constant input power, then the heat dissipation will be better.

Reference
[1] Liu S Luo X B 2011 LED Packaging for Lighting Applications: Design, Manufacturing and Testing Singapore Wiley
[2] Krames M R Shchekin O B Mueller-Mach R Mueller G O Zhou L Harbers G Craford M G 2007 J. Disp. Technol. 3 160
[3] Schubert E F Kim J K 2005 Science 308 1274
[4] Ma Y Hu R Yu X Shu W Luo X 2017 Int. J. Heat Mass Transfer 106 1
[5] Abdelmlek K B Araoud Z Ghnay R Abderrazak K Charrada K Zissis G 2016 Appl. Therm. Eng. 102 251
[6] Chhajed S Xi Y Li Y L Gessmann T Schubert E F 2005 J. Appl. Phys. 97 054506
[7] Ye H Mihailovic M Wong C K Y van Zeijla H W Gielen A W J Zhang G Q Sarroa P M 2013 Appl. Therm. Eng. 52 353
[8] Xie B Hu R Luo X 2016 J. Electron. Packag. 138 020803
[9] Nakamura S Fasol G 1997 Blue Laser Diode: GaN Based Light Emitters Lasers Berlin Springer-Verlag 16 20
[10] Xu H Y Chen X F Peng Y Xu M S Shen Y Hu X B Xu X G 2015 Chin. Phys. 24 067305
[11] Qi W J Zhang M Pan S Wang X L Zhang J L Jiang F Y 2016 Acta Phys. Sin. 65 077801 in Chinese
[12] Norazlina M S Dheepan Chakravarthii M K Shanmugan S Mutharasu D Shahrom M 2017 Chin. Phys. 26 098901
[13] Chen H T Y J Gao Y L Chen Z Zhuang R R Zhou X F Zhou H 2012 Acta Phys. Sin. 61 167104 in Chinese
[14] Arik M Becker C Weaver S Petroski J 2004 Proc. SPIE 5187 64
[15] Ma L Yang Y Liu J 2013 Heat Mass Transfer 49 85
[16] Tian L X Wen S S Yao R H Chen Y C Xie J N 2014 Acta Opt. Sin. 34 1123002 in Chinese
[17] Christensen A Graham S 2009 Appl. Therm. Eng. 29 364
[18] Joo Y H Kim S J 2015 Int. J. Heat Mass Transfer 83 345
[19] Ha M S 2009 Thermal Analysis of High Power LED Arrays MS Thesis Atlanta Georgia Institute of Technology